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ABSTRACT 

The radio emission from the youngest known Planetary nebula, SAO 244567, 
has been mapped at 1384, 2368, 4800, 8640, 16832 and 18752 MHz by using the 
Australian Telescope Compact Array ( ATCA) . These observations constitute 
the first detailed radio study of this very interesting object, as they allow 
us to obtain the overall radio morphology of the source and to compute, for 
the first time, the radio spectrum up to millimetre range. Radio emission is 
consistent with free-free from a wind-like shell, which is also the region where 
most of the [OIII] comes from as revealed by HST images. Physical parameters 
of the radio nebula and of the central star were derived, all consistent with 
SAO 244567 being a very young Planetary Nebula still embedded in the dusty 
remnant of the AGB phase. The optically thin radio flux density appear to 
decrease when compared to data from the literature. Even very appealing, 
the variability of the radio emission, probably related to the evolution of the 
central object, needs further investigations. 
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1 INTRODUCTION 



The fate of a star with a main sequence mass in the range from 1 to 8 solar mass is well 
established: it goes through the asymptotic Giant Branch (AGB) phase, then into the Plan- 
etary Nebula phase (PN hereafter) and eventually it will finish its evolution as a White 
Dwarf. However, the details of these evolutionary phases, such as the formation and the 
early evolution of PNe is far to be completely understood. It is not clear, for example, which 
kind of mechanism changes the more or less regular and symmetric circumstellar envelopes 
(CSEs) around PNe progenitors (AGB stars) into the quite complex morphologies observed 
in young and more evolved PNe. 

New clues on the process of PNe formation can be provided by the analysis of the physical 
characteristics of objects in the short phase between the end of the AGB and the onset of the 
ionization in the nebula. For this purpose, many authors have tried to identify very young 
PNe or proto Planetary Nebula (PPNe) but this has revealed to be quite difficult as this 
evolutionary phase is very rapid and because the central object is often heavily obscured by 
the thick circumstellar envelope formed during the AGB phase. 

Among post-AGB objects, SAO 244567 appears to be unique, as it is evolved so rapidly 
that strong spectral and total lu minosity changes were followed in a human life time scale. 
Parthasarathy fc Pottashl (119891 ) classified SAO 244567 as a post-AGB star on the basis of 
its high galactic latitude and because of its far-infrared (IRAS) colours similar to those of 
known PNe. Very soon it was realized that, in few decades, its optical and ultraviolet spectra 
have developed characteristics typical of the presence of a nebu la making SAO 244567 one 
of the youngest PN never discovered (jParthasarathy et al.lll993l ). 

In particular, its optical spectrum has evolved quite rapidly: reported by iHenizd (119761 ) 
with only H a in emission, more recentl y it has shown strong forbi dden nebular emission lines, 
which are consistent with a young PN (jParthasarathy et al.lll995l ). These results confirm that 
SAO 244567 has turned into a PN within the last 20 years and makes the source a perfect 
target for studying the early structure and evolution of Planetary Nebulae. 



SAO 244567 has been o 



Telescope (jBobrowsky et al 



jserv ed in 1996 with WFPC2 on board of the Hubble Space 



19981 ). These observations pointed out that most of the nebular 



emission originates from a ring/ellipse, whose major axes extends for ~ 1.6", and the presence 
of low density collimated outflows. The authors interpreted the ring structure as the remnant 
of the mass-loss during the AGB phase, while the low-density bipolar structures as results 



ATCA observations of the very young Planetary Nebula SAO 244567 3 



Table 1. Summary of results: flux densities are determined by fitting of visibility curves. 



Frequency 


Flux 


rms 


phase coX Flux 


MHz 


mjy 


mjy 


mjy 


Epoch 2000 








4800 


57.6 ± 1.7 


0.3 


4500 ± 50 


8640 


52.0 ± 1.6 


0.4 


3400 ± 50 


Epoch 2002 








1384 


36.6 ± 1.1 


0.2 




2368 


46.9 ± 1.8 


1.2 




4800 


48.8 ± 1.5 


0.2 


5234 ± 20 


8640 


46.6 ± 1.4 


0.2 


4271 ± 20 


16832 


43.8 ±2.0 


0.3 




18752 


42.8 ±2.0 


0.3 





from fast-wind experienced by the star during the post-AGB phase. The existence of a 
companion star, detected at 0.4 arcsec from the central star, indicates that the ring/eclipse 
is probably a circumbinary envelope and that binarity may play an important role in the 
shaping of PNe. These images are one of the best examples where nebular structures, that 
appear to collimate fast outflows, are evident. 

As typical for very young Planetary Nebulae (YPN), SAO 244567 shows a strong infrared 
excess and it is associated to the IRAS source IRAS 17119 — 5926. Measured non-corrected 
IRAS fluxes are 0.65, 15.50, 8.20 and 3.52 Jy at 12, 25, 60 and 100 /mi respectively 



In spite of the nume rous optical studies o: 



on its radio properties. iParthasarathy et al. 



this in teresting; object, very little is still known 



19931 ) briefly reported on ATCA 6 and 3 cm 



observations obtained in 1991. The measured flux densities are 63.3 ± 1.8 mJy and 51 ± 12 
mJy at 6 and 3 cm respectively. 

In this paper we present new ATCA observations of SAO 2444567 aimed to determine 
the radio properties of the nebula associated to this very young PN. 



2 OBSERVATIONS 

We carried out radio observations of SAO 244567 in March 2000 and August 2002. In both 
epoch the target was observed with the Australia Telescope Compact Array (ATCA) at 
Narrabri. The ATCA consists of six 22-m antenna, five of which lie along a three kilometer 
railway track, oriented east- west, and the other antenna lies on a similar railway track three 
kilometers away thus providing a 6 km maximum baseline lengths. 

As each ATCA antenna is equipped with dual-frequency systems, it is possible to observe 
simultaneously at two different frequency bands. 
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Figure 1. T he radio map of SAO 2 44567 obtained with the ATCA at 8.4 GHz (contour levels) superimposed with the HST 
image in H a l lBobrowskv et al]|l998h . The synthetic beam is shown in the lower left corner. 



2.1 The 2000 data 



In the first epoch, March 2000, the observations were simultaneously performed at 6 (4800 
MHz) and 3 cm (8640 MHz), with a 128 MHz bandwidth divided in 32 spectral channels, 
with an integration time of 30 sec. The observations were carried out in the 6D configuration, 
which, with a full synthesis, usually provides a typical beam size {0 syn ) of about 2" at 6 cm. 

The flux scale was fixed by daily observations of the primary flux calibrator 1934 — 638, 
whose assumed flux densities are 5.83 and 2.84 Jy at 4800 and 8640 MHz respectively. In 
order to achieve a good phase calibration, the phase calibrator 1718 — 649, which is 5.5 
degree away from SAO 244567, was observed at the beginning and at the end of each source 
observations. 
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At this epoch the source was observed in the contest of a survey program aimed to 
detect radio emission from a sample of post-AGB stars. Therefore, for a good compromise 
between UV coverage and total integration time, the observations were performed in snap- 
mode, which consists of short on source runs carried out at different hour angles. A typical 
observation of our target consisted of 7/8 cuts, each of 15 minutes long. In the case of SAO 
244567 we observed the source with 14 cuts, for a total of 3.5 hours. The observations were, 
however, obtained in different days starting from March 19 to March 22. 

Due to serious hardware problems at one of the antennas, the observations suffered of a 
degradation of the available angular resolution, resulting in a beam of 7.6" x 2.7" at 6 cm 
definitively worst than the angular resolution achievable with the 6D configuration. 



2.2 The 2002 data 

We have successively re-observed the source this time with a proper full tracking for better 
UV coverage. The observations were carried out on August 24 and August 29, 2002. In 
both dates the array was in 6C configuration. On August 24 we observe simultaneously at 
16832 and 18752 MHz, with a 128 MHz bandwidth divided in 32 spectral channels, using an 
integration time of 30 sec. At that time only 3 antenna were equipped with mm backends. 
High frequency observations are strongly affected by atmospheric phase fluctuations. To 
minimize this effect, we quickly switched between target and phase calibrator, with typical 
duty cycle of 6 min (four min on target and 2 min on phase calibrator) for a total of on 
source time of 7 hours. Reference pointing on our phase calibrator (1718 — 649) was also 
performed every hour of observing time. Flux scale was fixed by several observations of 
1924 — 638 and the flux of 1718 — 649 was boot-strapped using only those scan of primary 
calibrator performed at the same elevation as 1718 — 649. 

On August 29, the observations were carried out at 4800/8640 MHz and at 1384/2368 
MHz. The observations were performed following the same procedure as on August 24, using 
1718 — 649 as phase calibrator, with a typical on source scan of 15 minutes. However, only 
30% of the total on source time of 10 hours was spent observing at the lowest frequencies. 
The flux scale was fixed relatively to 1934 — 638 whose flux is assumed to be 14.94, 11.59, 
5.83 and 2.84 Jy at 1384, 2368, 4800 and 8640 MHz. 
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3 RESULTS 

The data from both epochs were reduced with the MIRIAD package, following the standard 
reduction steps. Visibilities were weighted by applying the robust weighting (robust param- 
eter^. 5). This methods allows to obtain the same sensitivity as with natural weight but 
with a much better beam shape. Finally the obtained maps were cleaned down to few times 
the theoretical noise, estimated to be of the order of 0.1-0.2 mJy 

The map with best combination of resolution and sensitivity is that obtained at 8640 MHz 



in 2002. The radio map (Fig. [T]) reveals a slightly exte nded structure, w 



rose overall size 



19981 ). However, 



(~ 1.7") compares quite well with that observed in Ha ( iBobrowsky et al. 
the limited ATCA angular resolution prevented us to point out any fine detail of the radio 
structure. 

Total flux density and the angular dimension can be usually estimated directly from the 
obtained map, by a 2D Gaussian fit of the radio source. However, because an interferometer 
will spatially filters out structures larger than the resolution corresponding to the minimum 
antenna spacing, the evaluation of flux density by Gaussian fit becomes less precise for 
partially resolved sources. This effect will be more evident at higher frequencies. In this 
cases it is advisable to derive flux densities directly from UV data, by fitting the visibilities 
at zero baseline. When the sources are angularly small, different methods to derive the flux 
density, i.e. from the map or from the fit of visibility data, provide consistent flux values to 
within a few %. In our case, differences between total fluxes as measured from the map and 
from fit at zero baseline are of the order of 0.1 percent up to 4800 MHz, but become more 
important at 18750 MHz where it goes up to ~ 17%. For this reason, in the following analysis 
we will considered the flux density derived from visibility fits. Results are summarized in 
Table I, where the observing frequency, the measured flux density, with its associated a 
and the rms of the visibility fit are reported. An error in the flux calibration can result in a 
systematic error in flux of the observed source. Therefore the a associated to the flux density 
estimation is derived from: 



a 



where rms is the error associated to the fit, cr ca i is the error associated to the absolute flux 
density scale, which is accurate to within few (3-5) percents and S is the derived flux density. 
In order to successively compare results from different epochs, for frequencies with multi- 
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Table 2. Parameters of different models for the radio nebula around SAO 244567. 
Parameter Value 



Model 1 (Wind Shell) 



Density 

Shell internal radius ijj nt 
Shell external radius -R cx t 
Shell Temperature 

Model 2 (Shallow Shell) 



n(-Rint) = 2.5 X 10 4 cm- 3 , n oc R~ 2 
0.75 arcsec 
1.5 arcsec 
9600 K 



Density 

Shell internal radius 
Shell external radius 
Shell Temperature 

Model 3 (Sphere) 



n = 1.45 X 10 4 cm" 3 
0.65 arcsec 
1.3 arcsec 
10000 K 



Density 

Sphere size 

Sphere Temperature 



n = 1.23 X 10 4 



1.4 arcsec 
10000 K 



epoch data we also report, in the last column of Table I, the flux density of the phase 
calibrator (1718 -649). 

We note that the flux density at 4800 and 8640 MHz appears to decrease from 1991 
(jParthasarathy et al.lll993l ) to 2002 (this paper). This decrement appears to be real and not 
related to a variation in the absolute flux scale as the flux density of the phase calibrator 
changes between the last two epochs but with a different trend, and it is actually increasing 
maintaining the same spectral index. In the following analysis only the multi-frequency 
epoch 2 data set will be used. 



4 THE NEBULA 

4.1 The Radio Morphology 

Some information on source morphology can be derived from the analysis of the visibility 
data, i.e. the fringe amplitude as function of the interferometer spacing. Since the ATCA 
is a linear interferometer, each of the different 15 minutes on source scans (cuts) will give 
a maximum resolution in one direction and the corresponding visibility will be function 
of the source morphology in that direction. Different visibility corresponding to different 
cuts look quite similar and this indicates that the overall morphology of the radio source 
is symmetric with respect to a central point. The small elongation in the S-W direction, 
probably reflects the UV coverage, and therefore the synthetic beam with which the visibility 
data are convolved and it is not intrinsic to the source. 
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We have modelled the visibilities of the entire observing run, averaged over 15 minutes, 
assuming 3 different symmetric morphology for the radio source, i.e.: 1) a shallow shell, 
0.75" thick, with a density decreasing with distance from the central star (n oc ^); 2) a 
constant density shallow shell, 0.65" thick; 3) a constant density sphere, with a radius of 
1.4". 

The best-fit parameters of the 3 models are summarized in Table II and results, relative 
to visibilities observed at 8640 MHz, are shown in Fig. [2j In the first panel the expected 
visibility from a stellar wind is shown (dotted line) for comparison. From this analysis we 
may conclude that the source morphology is consistent with both a shallow shell (with 
constant or with oc \ density) and with a constant density sphere. 



4.2 The physical properties 

Information on several parameters of a radio source can be derived by analyzing its radio 
spectrum, in particular when the transition region between optically thick and optically thin 
regime can be pointed out. The analysis of visibility data indicates a symmetric morphology, 
consistent with both a constant density sphere and a wind or constant density shell. Very 
little studies on modelling of free-free radio spectrum of PNe have been conducted, mostly 
because of the very few multi-frequency spectra available in the literature. Among them, 



Aaquist fc Kwokl (119911 ). by analyzing a sample of compact (young) PNe, noted a trend of 



increasing value of spectral index a (flux oc u a ), evaluated in the optically thick part of the 
spectrum, as the turnover frequency (v c ) decreases. This was interpreted as evidence, in the 
very early stages of nebula evolution, for most of the free-free emission coming from the AGB 
progenitor's wind, with typical a = 0.6. Taylor et al. (1987) modelled the radio continuum 
of 18 compact PNe with a shell model, with a radially dependent density. Observed radio 
spectra are well represented by a wind-shell, supporting the idea that free-free originates 
from a region formed with a constant stellar mass-loss during the precedent red giant phase. 

These precedent results together with the fact that SAO 244567 is a very young PN, will 
lead us to assume for SAO 244567 a wind-shell morphology. Model 1 well reproduces the 
observed visibilities at each frequency (Fig. [3]). The foreseen radio spectrum, superimposed 
to the observed data, is shown in Fig. HI with a turnover frequency, i.e. the frequency at 
which the nebula becomes transparent, between 2368 and 4800 MHz. The model will allow 
us to get an estimate of the total ionized mass contained in nebular shell. In the hypothesis 
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T 




Baseline (k A) 

Figure 2. Comparison between different modelling of the radio source. Both a shallow shell (with density oc 4j- or constant) 
and a constant density sphere are consistent with the observed visibilities. A stellar wind extending to infinity is shown as 
example. 
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Figure 3. The modelling of visibilities at various frequencies. Fits are obtained assuming a shallow shell, with decreasing 
density (oc r -2 ), with T = 9600 K and n = 2.5 X 10 4 cm 3 at the inner radius 



of a pure hydrogen nebula at 5.6 Kpc (jParthasarathy et al.lll993l ). we derive a total ionized 
mass of ~ 0.065 Mq. This value depends on the value of volum e filling factor of the nebula , 
which has been assumed to be 1 . A similar value was derived by iParthasarathy et al.l (119931 ) , 



while 



Bobrowsky! (119941 ) derived a value of 0.2Mq from the Hp flux. 



In the spectral region where the source is optically thin we can use the observed ra- 



dio flux density p 



( iTerzian fc Dickey 



us th e angular size to derive the mean emission measure as follows 



1973|): 



< EM >- 



J u EMdVt 5.3 x 10 5 F 8 . 4 G h z 



1.1 x 10 



7 



n e 2 

where F 8 . 4 qh z is the measured radio flux density, in mJy, at 8.4 GHz, and 9, in arcsec, 
the angular dimension of the radio emitting region as derived from the model. The mean 
emission measure is expresses, as usual, in cm _6 pc. 



Young PNe should have emission measures of the order of 10 6 -10 8 cm 6 pc ( ITerzian fc Dickey 



1973 



Kwok et al. 



19811). 
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Figure 4. The radio spectrum of SAO 244567 obtained from the fit showed in figure taking as flux density the zero spacing 
value. 

From the radio optically thin flux we can also derive the excitation parameter (U exc ) 
required to account for the measured radio flux: 

(1) 

where F v is the optically thin radio flux density, expressed in Jy, at the observing frequency 
v (in GHz); T is the nebula temperature, expressed in 10 4 K, and D-^ pc is the source distance 
in Kpc. 

As the PN associated to SAO 244567 is very young we can assume that it is ionization 
bounded; in this case the excitation parameter is directly related to the number of ionizing 
photons emitted by the central star. 

L uv = 1.23 x 10 5f W c 3 xc photons s" 1 (2) 

where (3 is the hydrogen recombination co efficient summ ed over all levels above the ground 
level (~3x l(T 13 cm 3 s^ 1 for T ~ 10 4 K) ls P itzer!ll978k 

From equation (TjQ) and equation (j2J), assuming a nebular temperature of T ~ 9600 K and a 
distance of 5.6 Kpc, from the radio flux density measured at 4.8 GHz, we derive a number 
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of ionizing photons of L uv ~ 1.5 x 10 47 (photons s _1 ). 

This corresponds to the flux of Lyman con tinuum photon s of a BO-BO. 5 V with an effective 
temperature of T cS ~ (2.8 ± 0.2) x 10 4 K Jpanarialll973k 

This value for the effective temper ature of the central s t ar is s omehow lower than the value 
of T e g ~ 3.7 x 10 4 K estimated by 



Parthasarathy et al 



(119931 ) from the analysis of the UV 



spectra and corresponding to a 08V star. This may indicate the presence of significant 
quantity of dust well mixed with the gas in the inner regions of the nebula. Dust, in fact, 
competes with the gas for the absorption of the ionizing photons and for a given number of 
ionizing photons emitted by the central star, it will reduce the observed radio flux. 



THE DUSTY ENVELOPE 



One of the most important results from PNe studies based on IRAS data was the realization 
that, on the average, in young PNe about 40% of the emergent flux is emitted in the far- 
infrared (jZhang fc Kwoklll99ll ). This is due to the presence of an extended dusty envelope, 
around YPNe, which is the remnant of the precursor's wind not yet dispersed. In YPNe 
the spectrum emitted by the dust, bet ween 10 and 100 ton, i s usually well described b y 



a single temperature black body curve (jZhang fc Kwokl Il99d ; iStasinska fc Szcerbal Il999l ). 
We therefore derive a dust temperature of T dust = 137 ± 2 by using a least-squared fitting 
procedure to th IRAS measurements. This procedure provides, as a by-product, the total 
far- infrared flux (Fir), obtained by integrating, over the IRAS band (25 to 100 /zm), the 
Planck curve that fits the IRAS data. 

It is therefore possible to define the far infrared excess (IRE) as the ratio of the observed 
total far infrared flux (Fir) over the expected total infrared flux. Under the hypothesis that 
the far infrared flux is due to thermal emission from dusty grains heated by Lya photons, this 
ratio is unity. However, in y oung PNe heati ng by direct starlight is important and IRE can 
be much higher than unity. iPottaschl (119841 ) has derived a formula to express the expected 
total infrared flux in terms of optically thin radio flux density: 
Fm 



IRE = 1.07- 



2.8 



-^8.6 GHz 

where Fir is expressed in 10 -14 Wm -2 , -Fs.eGHz in niJy, and the high density approximation 
has been adopted. This value is consistent with the hypothesis that SAO 244567 is a very 
YPN where the dust plays an important role. 
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To qu antify the total mass of the dust we will use the optically thin expression (IHildebrand 



1983|): 



M, 



dust 



FQ0fj, m D 2 



xBv(T dust ) 

where F 60Mm is the flux density measured at 60//m, D is the distance to the source (Kpc), 
Bv(Td nst ) is the value of the BB function of Td US t, evaluated at 60//m and x is the mean 
absorption coefficient for the grains (cm 2 g -1 ). This formula provides only an estimation of 
the total dust mass as a single temperature and an unique size for the grains is assumed. 
Moreover, while this dust mass determination is independent on the distribution of the 
grains, it strongly depends on the chemical composition of the grains. 
There are still problems in understanding the chemical composition of the dust component 
in PNe. Usually, in carbon rich PN, dust is mainly composed of carbon-based grains, while 
in O-rich PN, dust is mainly composed of different forms of silicates. As for SAO 244567 we 
don't have any information on its chemical composition, we compute the dust mass for two 
different value of x, namely 53.45 cm 2 g _1 , whic h corresponds to circumste llar silicates and 
145.32 cm 2 g _1 , which corresponds to graphite (jStasinska fc Szcerballl999l ). This results in 
a total mass dust of 2 x 1O _4 M , in case of silicates, and 7.5 x 1O _5 M , in case of graphite. 
If we assume that dust and gas occupy the same volume, we may derive a dust to gas ratio 
of 3 x 10~ 3 and 1.15 x 10~ 3 for silicates and graphite r espectively. These values a re in the 
range for M dust M gas derived for a large sample of PNe ( Stasinska fc Szcerba Il999 ). 



6 THE VARIABILITY OF RADIO EMISSION 

One of the most striking and unexpected characteristic of radio emission from SAO 244567 is 



i ts var iability. In Fig.[5l we report the flux density at 4800 MHz, as measured by lParthasarathy et al 



( 119931 ) and in this paper, indicating a linear decreasing trend of the optically thin flux. The 
dashed line is the weighted fit to the 4800 MHz data points, which results in a flux decrement 
of ~ 1.3 mJy year -1 Assuming that the radio emission from the source will keep this trend a 
flux density of ~ 43 mJy in 2008 is foreseen. In the same figure we also include the 2000 and 



2002 8 640 MHz data, without considering the quite noisy datum from 



(1993) 



Parthasarathy et al 



In the hypothesis that, between March 2000 and August 2002 (AT), the ionized mass re- 
mains the same (i.e. any effects due to the evolution of the central object are not considered), 
and the nebula is expanding homologous, we have: 
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M Tot oc 

^2000-^2000 ~~ n 2002-R2002 



(3) 



Starting from the results for epoch 2002, from equation ([3]) we obtain constraints for 
epoch 2000. To reproduce the flux density observed at 4800 and 8640 MHz in epoch 2000, 
we then used the model, as illustrated in par 4.1, obtaining a reasonable fit that satisfies 
equation ([3]), for T e s = 9600 K, as in 2002 fit, with an internal radius of 0.71", an external 
radius of 1.4" and a density, at the base of the shell, of n2ooo — 2.9 10 4 cm" 3 . To obtain such 
an expansion we need the nebula to expand in a time AT with an average velocity of ~ 850 
km sec -1 . We should stress here that a variation of radius of the order of few 0.1" would be 
not detectable with present instrumentation (ATCA). 

Hig h expansion veloci ty have been observed is the ultraviolet lines of some PPNe and 
YPNe (jGauba et al.ll200ll ) and are usually related to fas t -wind s events occurring during the 
post-AGB phase. For SAO 244567 iParthasarathy et al.l (119951 ) reported, in 1988, an ultra- 
violet spectrum characterized by strong P-Cyg profiles, indicating the presence of a strong 
wind with expansion velocities of the order of ~ 3000 km sec -1 . This spectral structures 
were observed until 1993 and by 1994 the y were almost vanish e d, ind icating the end of a 
rapid mass- loss event. At the same times, IParthasarathy et al.l (119931 ) derived, from high- 
resolution [OIII] profile, an expansion velocity of the order of ~ 8 km sec -1 , consistent with 
the fact that the ionized shell, from which most of the [OIII] comes from, is the remnant of 
mass-loss during the previous AGB phase. Both results are consistent if we assume that UV 
lines are produced in the low-density polar region, where episodic rapid mass-loss event can 
occur, while the bulk of the ionized material, and also the radio emission, is localized in the 
main ring/shell, as indicate by the HST images, which is expanding with typical velocity 
as expected for a AGB wind. Therefore, we may conclude that an expansion of the main 
shell can not explain the observed decrement in the radio flux density because the required 
expansion velocity is too high. 

High velocity mass-loss events may, however, contribute to th e total ionization of the 



nebula. This is what is claimed by 



Sanchez Contreras et al. 



(120041 ) to explain the decrease 



in millimetic, free-free optically thin flux observed in the PPNe CRL 618. According to these 
authors, in this source there is a substantial contribution of shocks, originating from inter- 
action between the fast post-AGB wind with the CSE environment, to the total ionization. 
As the fast wind stops, the ionization rely only on the central star radiation field and the 
optically thin radio emission is, therefore, decreasing. However, this explanation requires 
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episodic high-velocit y mass-loss events to justify the flux increasing reported by the same 



authors 2 year later (ISanchez Contreras et al 



sity has been also noted in NGC 7027 (jZijlstra et al 



2004b ). A decrease in the optically flux den- 



20071 ) and it is interpreted as related to 



the evolution of the central object. NGC 7027 is located just at the tip of the WD cooling 
track, where its luminosity starts to decrease while its temperature is still increasing. This 
leads to a decreasing number of ionizing photons. Since the optically thin free-free emission 
is directly related to the ionizing flux from the central object, this would imply a decreasing 
optically thin radio flux density. 

SAO 244567 has shown many spectral changes in the last 20 years. Most notably, the changes 
in the UV flux level, which suffered a factor 2.83 decre ment in 7 years, while the central was 
gradually becoming hotter (jParthasarathy et al.lll995l ). implying a total drop of the central 
star luminosity. This pl ace the central star of SAO 244567 in the same position on the HR 

200 71), just on the top of the cooling WD track. How- 



diagram as NGC 7027 flZiilstra et al 



ever, as already pointed out by 



Parthasarathy et al.l (119951 ) ; iBobrowsky et al.l (119981 ) , such 



a very fast evolution is very difficult to be explained in the fra mework of actual post-AGB 



evolutionary models, that for a core mass object of ~ 0.056M Q ( jParthasarathy et al. 
foresee a much slower evolution. 



1993 ). 



7 SUMMARY 

We have presented ATCA multi-epoch, multi-frequency observations of the YPNe SAO 244567, 
aimed to derive its radio characteristics to complete the picture of this quite intriguing ob- 
ject. The best radio map, for sensitivity and resolution, has been obtained at 86400 MHz (3.6 
cm), which reveals a slight extended radio structure. The angular resolution of ATCA does 
not allow to evidence the fine details of the nebula as shown by HST observations. However, 
some morphological information have been derived by a model to the observed visibilities 
data, which are consistent with a wind-like shell whose external radius is in agreement with 
that of the main structure observed by HST. 

The mean emission measure and the infrared excess, as derived by our observations, are 
consistent with a very young planetary nebula, still embedded in its dusty envelope, remnant 
of the earlier AGB phase. By comparison between total ionized gas mass, as derived from 
the radio, and the total dust mass, as derived from the IRAS data, a dust to gas ratio has 
been derived, by considering two different types of chemistry for the dusty envelope. 



16 G. Umana, C. Trigilio, L. Cerrigone, C.S. Buemi and P. Leto 







70 



60 



50 



40 



30 




SAO 244567 

a = -1.27 +/-0.2 



1990 



1995 



2000 

Year 



2005 



2010 



Figure 5. The behaviour of flux density of SAO 244567 as measured at 4800 MHz, filled dots. Data are from lPart hasarathv et all 
jl993h and this paper. A weighted fit to the data is also shown, which foresees a radio flux density at this frequency of 43 mjy 
at the beginning of 2008. Measures at 8640 MHz are also plotted, diamonds. 



When compared with previous observations, the radio flux appears to vary. In particular, 
the decrement of the 6 cm flux (from 1991 to 2002), does not agree with an expansion of 
the main shell but appear to be related to the characteristics of the stellar object, whose 
fast evolution is difficult to be explained in the framework of the recent models of post-AGB 
evolution. Further multi-frequency monitoring and mapping with the ATCA fully equipped 
with millimetre receivers are necessary to confirm such variation with a longer time baseline. 
Up to now, variation (decrement) of the optically thin radio emission have been observed 
only in other two objects: the proto-PN CRL 618 and the PN NGC 7027. In both cases, the 
observed variations have been interpreted in terms of evolution of the central object. The 
paucity of objects where such variations have been followed makes further radio observations 
of SAO 244567 very important, as typical timescale of such variations can be used to test 
current evolutionary models. 



ATCA observations of the very young Planetary Nebula SAO 244567 17 
REFERENCES 

Aaquist O.B., Kwok S., 1991, ApJ, 378, 599 
Bobrowsky M., 1994, ApJ, 426, L47 

Bobrowsky M., Sahu K.C., Parthasarathy M., Garca-Lario P., 1998, Nature, 392, 469 

Gauba G., Parthasarathy M., Nakada Y., Fujii T., 2001, A&A, 373, 572 

Henize K.G., 1976, ApJS, 30, 491 

Hildebrand R.H., 1983, QJRAS, 24, 267 

Kwok S., Purton C.R., Keenan D.W, 1981, ApJ, 250, 232 

Panagia N., 1973, AJ, 78, 929 

Parthasarathy M., Pottasch S.R., 1989, A&A, 225, 521 

Parthasarathy M., Garcia-Lario P., Pottasch S.R., Manchado A., Clavel J., de Martino D., 
van de Steene G.C.M., Sahu K.C., 1993, A&A, 267, L19 

Parthasarathy M., Garcia-Lario P., de Martino D., Pottasch S.R., Kilkenny D., Martinez 
P., Sahu K.C., Reddy B.E., Sewell B.T., 1995, A&A, 300, L25 

Pottasch S.R., 1984, ASSL, 107, Planetary nebulae - A study of late stages of stellar evo- 
lution. Dordrecht, D. Reidel Publishing Co. 

Sanchez Contreras C, Bujarrabal V., Castro-Carrizo A., Alcolea J., Sargent A, 2004, ApJ, 
617, 1142 

Sanchez Contreras C, Sahai R., 2004, ApJ, 602, 960 

Spitzer L., 1978, Physical Processes in the Interstellar Medium. Wiley-Interscience, New 
York 

Stasinska G., Szczerba R., 1999, A&A, 352, 297 

Taylor A.R., Pottasch S.R., Zhang C.Y, 1987, A&A, 171, 178 

Terzian Y., Dickey J., 1973, AJ, 78, 875 

Zhang C.Y., Kwok S, 1990, A&A, 237, 479 

Zhang C.Y., Kwok S, 1991, A&A, 250, 179 

Zijlstra A., Per ley R., van Hoof P., 2007, NRAO Newletters, 111, 7 



